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Abstract 

We study composite systems of light {p and uj) and heavy (D*) vector mesons 
^SJ ! by using the interaction within the hidden gauge formahsm. We find a strong at- 

^ I traction in the isospin, spin channels {I,S) = (1/2,0); (1/2, 1); (1/2,2) with positive 

■ parity. The attraction is large enough to strongly bind these mesons in states with 
, these quantum numbers, leading to states which can be identified with D2(2460) and 

Q"^ I probably with I?* (2640), the last one without experimental spin and parity assign- 

■ ment. In the case of I = 3/2 one obtains repulsion and thus, no exotic mesons in 
this sector are generated in the approach. 

• ^ - 

X 

5^ ; 1 Introduction 

Recently a study of the pp interaction with the hidden gauge formalism [HE], [3] was carried 
out in [3]. The hidden gauge symmetry (HGS) was introduced by Bando-Kugo-Yamawaki 
where the p meson was regarded as a dynamical gauge boson of the HGS of the non- 
linear sigma model. One of interesting facts which is relevant in the present discussion 
is that there is a strong attraction in the isospin and spin channels (/, S) = (0, 0) and 
(0, 2), which is enough to bind the pp system leading to a tensor and a scalar meson which 
could be identified with the /o(1370) and /2(1270) meson states [4j. In a later work the 
radiative decay of these states in 77, a channel with rates very sensitive to the nature of the 
resonances, was studied [5] obtaining results in agreement with the PDG |^ for the case 
of the tensor state and in qualitative agreement with preliminary results at Belle for the 
scalar state [7j. The work of [3] has been recently extended to SU(3) for the interaction 
of the vectors of the nonet, where several states, which can be identified with existing 
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resonances are also dynamically generated [8]. The success of the approach encourages us 
to study the charm sector in the lightest case, the one for the interaction of the p, uj and D* 
mesons. Another possible approach to this work could be done following the lines of [9] for 
meson-baryon interaction involving SU (8) symmetry. Work in this direction is in progress 
[To] but we can advance that while for pseudoscalar-baryon interaction the approaches are 
equivalent, when vector meson are involved the results are very different [11]. 

The starting point in our approachis the interaction of vector mesons among them- 
selves provided by the hidden gauge formalism, which now has to be generalized to SU(4) 
to accommodate the charm vectors D* into the framework. Admitting that SU(4) is more 
strongly broken than SU(3), the SU(4) symmetry is invoked in the basic hidden gauge 
Lagrangians but is already broken in the vector exchange diagrams that provide the am- 
plitudes for the vector- vector interaction. What we find is a strong attraction in the isospin, 
spin channels {I,S) = (1/2,0); (1/2,1); (1/2,2); which leads to bound p{uj)D* states in 
all these channels. In the case of / = 3/2 we find repulsion and hence we do not generate 
states that would qualify as exotic from the qq picture. The states that we find qualify as 
mostly pD* molecules, and fit nicely with the experimental states 1^2(2460) and D*(2640). 
The present study would, thus, suggest for these states a different nature than the one 
usually assumed in terms of qq states [HI [121 113] . 



2 Formalism for VV interaction 
2.1 Lagrangian 

We follow the formalism of the hidden gauge symmetry (HGS) for vector mesons of [Tl[2] (see 
also [H] for a practical set of Feynman rules). The Lagrangian involving the interaction 
of vector mesons amongst themselves is given by 

cju = -\{v,.vn , (1) 

where the symbol () stands for the trace in the SU{4) space and V^jy is given by 

V^. = d^V, - d,V^ - tg[V^, K] , (2) 

with g given by 

9 = ^, (3) 

and / = 93 MeV the pion decay constant. The value of g of Eq. ([3]) is one of the ways 
to account for the KSFR relation [TJ] which is tied to vector meson dominance [TB]. The 
vector field is represented by the SU{A) matrix which is parameterized by 16 vector 
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Figure 1: Terms of the Cui Lagrangian: a) four vector contact term, Eq. (jSj); b) three- 
vector interaction, Eq. (jSD; c) t and u channels from vector exchange; d) s channel for 
vector exchange. 

mesons including 15-plet and singlet of S'f/(4), 

/ 4 + ^ p+ K*+ D*^ \ 



V2 V2 

p- -^ + ^ K*^ 

K*- D* 

V D*+ D*+ J/tp J 



(4) 



where the ideal mixing has been taken for u, and J/ip. The interaction of Cjjj gives rise 
to a contact term coming for [V^, V^] 

^ni = ^^{V.V.V^V' - V^V^V^V) , (5) 
depicted in Fig. [Ha), and on the other hand it gives rise to a three vector vertex 

4//^ = M{d,V, - duV.W^V) , (6) 

depicted in Fig. [T]b). This latter Lagrangian gives rise to a VV — > VV interaction by 
means of the exchange of one of the vectors, as shown in Fig. [T]c). 

The SU{A) structure of the Lagrangian allows us to take into account all the channels 
within SU{A) which couple to certain quantum numbers. In the present work we shall 
present results for the case of the pD* interaction. The formalism is the same used in 
[1]. Some approximations were made there which make the formalism handy and reliable, 
by neglecting the three-momentum of the vector mesons with respect to their masses. 
It is interesting to see that with this approximation one obtains [H] from the hidden 
gauge approach the chiral local Lagrangians which are used to study the interaction of 
pseudoscalar mesons among themselves and the pesudoscalar mesons with vector mesons 
and with baryons [T7[ [H] . 

2.2 Four- vector contact interaction 

Starting with the Lagrangian of Eq. ^ we can immediately obtain the amplitude of, for 
instance, p^D*^ p^D*^ corresponding to Fig. [2], in the particle base. 
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Figure 2: Contact term of the pp interaction. 

where the indices 1, 2, 3 and 4 correspond to the particles with momenta ki,k2, and k^ 
in Fig. [21 It is straightforward to write down all amplitudes for other channels. 

In the approximation of neglecting the three-momenta of the vector mesons, only the 
spatial components of the polarization vectors are nonvanishing and then one can obtain 
easily spin projection operators [4] into spin 0, 1, 2 states, which are given below: 

V^'^ = {i(e^6.e'^e'^ + 6,e.6^e'^)-le^eV}, (8) 

where the order 1, 2, 3, 4 in the polarization vectors is understood. We can then write the 
combination of polarization vectors appearing in Eq. ([7]) in terms of the spin combinations 
and thus we obtain the kernel of the interaction which will be later on used to solve the 
Bethe-Salpeter equation. However, it is practical to construct the isospin combinations 
before the spin projection is done. 

Recalling that we have an isospin doublet with {—D*^, D*^) and one isospin triplet, 
{p^,p^, — p^), the I = 1/2 and 3/2 combinations are written as 



|pD*,/=l/2,/3 = l/2) = ^^|p+D*°)-i=|/D*+), 

\pD*,I = 1/2, h = 3/2) = -L|p+D*°) + ^|p°Z^*+) . (9) 

We then find the amplitudes in the isospin base by forming linear combinations of the 
amplitudes in the particle base weighted by the Clebsh-Gordan coefficients as given in 
Eq. m, 
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These amplitudes, after projection into the spin channels, give rise to the following kernels 
(potential) for J = 1/2, 

^(7=1/2,5=0) ^ ^5^2^ 
(7=1/2,5=1) _ ,92 



and also for the case of / = 3/2, 



^(/=3/2,5=0) ^ _4^2^ 
^(/=3/2,5=l) ^ 

^(7=3/2,5=2) ^ ^2^7^ (12) 

When one or two p meson(s) are replaced by the u meson(s), where is only one isospin 
state / = 1/2, we find the following interaction terms. 



(7=1/2) _ a u _ u u _ u u-s 

1 



2 



After projection in spin they become for pD* uoD* 



(7=1/2,5=0) _ _ 2 
(7=1/2,5=1) _ ,3^3 2 



t^;^sjif = +^/. (14) 



and in the same way, we have for ujD* ujD*, 



(7=1/2,5=0) _ _ 2 

^^0*^1^0" — y 1 

(7=1/2,5=1) _ S o 

^a;D*-^a;7)* — ^ 



2.3 /o meson exchange terms 



From the Lagrangian of Eq. ^ we get the three- vector vertex as depicted in Fig. [31 For 
practical purposes it is convenient to rewrite the three-vector Lagrangian of Eq. ([6]) as, 

= ig{{V^d,V^ - d^V^V^)^) . (16) 
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D*^{k-q,e') 



D*\k,e) 



Figure 3: three- vector vertex diagram. 

In Eq. (fT6l) we have a three-vector vertex, where any of the three vector fields can corre- 
spond in principle to the exchanged vector in the diagram of Fig. 1(c). Nevertheless, by 
assuming that the three-momenta of the external vectors can be neglected as compared 
with the vector mass, the polarization vectors of the external vector mesons have only 
spatial components. Then by looking at the Lagrangian of Eq. (|T6l) we see that the field 

cannot correspond to an external vector meson. Indeed, if this were the case, the v 
index must be spatial and then the partial derivative is replaced by a three-momentum 
of the vector mesons which is neglected in the approach. Then corresponds to the 
exchanged vector and this simplifies the calculation. The approximation of neglecting the 
three-momenta of the external vectors corresponds to the consideration of only the s-wave. 

The vertex function corresponding to the diagram of Fig. [3] is given by 

+ (-^ V- + zfc,e^)e(°) V^} . (17) 

With this basic structure we can readily evaluate the amplitude of the first diagram of 
Fig. m to obtain 

i 

{h - hf - M2 + le 

-j^im - A:4).6(°) - z{h - hUf)e^^)^e^^> 
{-ik,,e^^^+ik,J^^)e^'^^e^'^^} . (18) 



- it 
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P" 





p\h - h) 




p' 




1*0 



Figure 4: Vector exchange diagrams for pD* pD*. 

Recalling that the three-momenta of the external particles is small and neglected, we arrive 
at the following expression: 

t^;:lo^,+ D'0 = {h + h) ■ {h + h) 6,e.6^6- . (19) 

By looking at the structure of the second diagram in Fig. H] we find the following result for 
the amplitude: 

2 

ti+U^pOn^+ = V2j-{h + h) ■ ih + h) e^e.e^e^ (20) 

p 

We note that the amplitude p^D* —>■ p^D* with p^ exchange vanishes because the three- 
p° vertex does not exist due to isospin invariance. We can see that in all the cases the 
combination of vector polarizations is the same. The isospin projections give us 

t%^Zl& = -^ih + h).ih + h)e,e,e'^e^ . (21) 

Now using the equations for the spin projections we can split the terms into their spin 
parts and we obtain 

t%'/Z%y^'^'''' = ^^{h + k,).ik, + k,), (22) 
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that is: we find spin degeneration in the amphtudes which involve the exchange of one p 
meson. These structures can be simphfied using momentum conservation and one finds: 



(ex,7=l/2,S=0,l,2) _ 9 ™2 2 (m^ m%,) 



P 

^pD.^pD* + }• (23) 

Note that for the exchange of a vector we do not have contribution from the ujD* channeL 
Indeed, the vertex uoju and ppu violate G-parity, and the the pujuj violates isospin. In 
the next section we consider the amplitudes which include also the exchange of one heavy 
vector meson, but we anticipate that the amplitudes calculated so far are more relevant 
than the other ones. 



2.4 D*- exchange terms 

In this section we are going to take into account the exchange of one heavy vector meson, 
D* or Z)*, in the channels pD* and uoD*, by means of the diagrams that we draw in Fig. [5l 
Note that in Fig. [5] the vertex uouoD* does not appear because it violates isospin. By means 
of the Lagrangian of Eq. ([6]) we arrive at the following amplitudes for the first and second 
diagrams depicted in that figure, 

D* 1 9^ 

tf+D*Tlp,D*+ = -^J^^ (ki + h) ■ {k2 + h) e^eye^" 

t%-.TlpOo^. = l^^{ki + h)-{k2 + h)e,eye^ . (24) 



The isospin decomposition can be done as before by making weighted sums of Eq. flMl) . 
We find, 

tfoCli^""^ = ^{k, + k,)-{k, + k,)e,eye^. (25) 
where k = M^/M^*. Upon spin projection we find, 

t%:-j;i^T'^'''-'''^ = -l'j^{k, + k,).{k, + k,) 

t%r:7^"'''^'' = l^ih + k,).ik, + k,), (26) 



8 



and similarly for I = 3/2, 



t 



(D*- ex, I =3/2,8=0,2) 



t 



{D*-ex,I=3/2,S=l) 



(/ci + k^) ■ {k2 + k^) . 



(27) 



Again by neglecting the three-momenta of the external particles and by taking the center- 
of-mass reference system, one can express the factors which involve momenta as 



{ki + /cg) ■ {k2 + ki) 
{ki + k^) ■ {k2 + ks) 



2^-^p 



m 



D* 



D* 



m 



D' 



+ 



2s 



2s 



(2J 



As we can observe, the spin degeneracy seen in the p-meson exchange amplitudes is lost 
in the D*-exchange amplitudes. However, it is broken only a little due to the suppressing 
factor K = m^/m^, ~ 0.15. 

The results are summarized in the columns titled "Z?*-exchange" in Tables [1], [2] and [31 
These new terms represent corrections of the order of 10% of the p-exchange ones. As can 
be observed in the total amplitudes, we find attraction in the sector I = 1/2, whereas the 
sector J = 3/2 turns out repulsive. It is interesting to see that the exotic 1 = 3/2 channel 
has a repulsive interaction. This seems to be a rather universal in this kind of studies 



I 


S 


Contact 


p-exchange 


Z)*-exchange 


~ Total[/(J^)] 


1/2 





+5g^ 


-2^ (A;i + fca) ■ {k^ + k^) 


-\l^{k, + ki)-{k2 + k,) 


-16^?2[l/2(0+)] 


1/2 


1 




-2^ {ki + k^) ■ {k2 + k^) 


+ ImJ (^1 + ^4) ■ {k2 + k^) 


-14.5(?2[l/2(l+)] 


1/2 


2 




-2^ (A;i + fcg) ■ {k2 + k^) 


-\^{k, + ki)-{k2 + k,) 


-23.5^?2[l/2(2+)] 


3/2 





-Ag^ 


+ £2 {ki + ks) ■ {k2 + ki) 


+ {ki + ki) ■ {k2 + fcg) 


+8<72[3/2(0+)] 


3/2 


1 





+ £2 (ki + ks) ■ {k2 + /C4) 


-i| (fcl + ki) ■ {k2 + ks) 


+8^7^[3/2(l+)] 


3/2 


2 


+2/ 


+ ^ (ki + k^) ■ {k2 + ki) 


{ki + ki) ■ {k2 + k^) 


+145^[3/2(2+)] 



Table 1: V{pD* pD*) for the different spin-isospin channels including the exchange of 
one heavy vector meson. The approximate Total is obtained at the threshold of pD*. 



2.5 T-matrix 

The results obtained in Tables [H, [2] and [3] provide the kernel or potential V to be used in 
the Bethe-Salpeter equation in its on-shell factorized form. 
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/ 


S 


Contact 


p-exchange 


D* -exchange 


~ Total [J (J^)] 


1/2 









+'^j^{h + k,)-{k2 + ks) 


0[l/2(0+)] 


1/2 


1 


I 3v^ 2 




-fj^{h + k,)-{k2 + ks) 




1/2 


2 






+fj^ih + k,)-{k2 + ks) 


^^2[l/2(2+)] 



Table 2: V{pD* ujD*) for the different spin-isospin channels including the exchange of 
one heavy vector meson. The approximate Total is obtained at the threshold of pD* . 



I 


S 


Contact 


p-exchange 


-exchange 


~ Total[/(J^)] 


1/2 











0[l/2(0+)] 


1/2 


1 






-li^ih + k,)-{k2 + ks) 




1/2 


2 


+ \9' 




+ \% (fcl + k^) ■ {k2 + k,) 


|^2[l/2(2+)] 



Table 3: V{u!D* ujD*) for the different spin-isospin channels including the exchange of 
one heavy vector meson. The approximate Total is obtained at the threshold of pD* . 



for each spin-isospin channel independently. Here G is the two meson loop function 

^_ ■ f d^Q 1 I /onx 

' J (27r)4 q^-ml + it [P - - ml + ie ' ^ ' 

which upon using dimensional regularization can be recast as 



^/s ^ —s + m2 — mi + Zpy's —s — 1712 + rrii + lj)^s' J 

where P is the total four-momentum of the two mesons, p is the three-momentum of the 
mesons in the center-of-mass frame. Analogously, using a cut off one obtains 

where Qmax stands for the cut off, Ui = iqi+mfy^'^ and the center-of-mass energy (P°)^ = s. 

Equation fl29l) inJ = l/2isa2x2 matrix equation with the amplitudes pD —>■ pD, 
pD —>■ ujD and ujD ujD in the elements (1, 1), (1, 2) and (2, 2). 

The formalism that we are using is also allowed for s-channel p or D exchange and we 
can have the diagram of Fig. El But we found in [1] that this leads to a p-wave interaction 
for equal masses of the vectors, and only to a minor component of s-wave in the case of 
different masses [8]. 
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Figure 6: S'-channel D* exchange diagram. 

3 Convolution due to the p mass distribution 

The strong attraction in the I = 1/2 and S = 0,1,2 channels will produce pD* bound 
states with no width within the present treatment so far. However, this is not strictly true 
because the p meson has a large width or equivalently a mass distribution that allows the 
states obtained to decay in pD* for the low mass components of the p mass distribution. 
To take this into account we follow the traditional method which is to convolute the G 
function for the mass distribution of the p meson [H] replacing the G function by G as 
follows 

1 rimp+2Tp)'^ _ ]^ I 

'^(^^ = m L ^dml{ )Im— G{s,ml,ml,) , (33) 



with 



N = / dml{—)Im— . , (34) 

J{mp-2Tp)^ Tc nil — rri'' + iL mi 



where Tp = 146.2 MeV and for F = F(m) we take the p width for the decay into the pions 
in p-wave 

F(m) = F,( ^ ; f'e{m - 2m.) . (35) 
mi — Ami 

p TV 

The use of G in Eq. fl29p provides a width to the states obtained as we will see in the next 
section. 



4 Results 

When one introduces the amplitudes obtained in Tables [H [2] and [3] as a kernel V in Eq. fl29l) . 
one finds bound states with zero width in the three different cases 1 = 1/2 and = 0, 1, 2. 
The pole positions are given in Table |H In Eq. fl3Tl) we have fixed the value of p as 
1500 MeV and we have fine tuned the subtraction constant, a, around its natural value of 
—2 [21J in order to get the position of the 5 = 2 resonance at its value of the PDG. To 
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Squared amplitude for 1=1/2 and S=0 
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Figure 7: 


Squared amplitude for I = 


1 /2 and 5* -- 


= 0, 



Squared amplitude for 1=1/2 and S=l 



2580 2600 2620 2640 
[MeV] 



2660 



p-mass distribution. 



quantify the freedom one has in this fine tuning we quote that the value of the mass that 
we obtain using a = —2, is 2346 MeV . The value of a for the pole positions given in the 
Table Sis -1.74. 



/ 


S 


v/i (MeV) 


1/2 





2592 


1/2 


1 


2611 


1/2 


2 


2450 



Table 4: Pole positions for the three different cases 

As it was explained in the previous section, the loop function G must be convoluted to 
take into account the width of the p meson. When we do it, we find bound states with a 
small width as we can see in Figs. [7] and [HI In fact, the widths found are ~ 5 MeV, 4 MeV, 
and MeV for S' = 0, 1 and 2, respectively. 

In Eq. ( !29l) . the amplitude close to a pole looks like 

T,, ^ , (36) 

Z- Zr 

where Re zr gives the mass of the resonance and Im zr the half width. The constants gi 
{i = pD* , ujD*), which provide the coupling of the resonance to one particular channel, 
can be calculated by means of the residues of the amplitudes as given in Table [51 In the 
PDG J6j, two states are listed, D*(2640) and D*(2460) in the sectors /(J^) = l/2(?-) and 
J(J^) = 1/2(2"'"), respectively. Comparing with the present model predictions as listed in 
Table [U we attempt to identify them with those states of = 1 and S = 2, respectively. 
The reason to identify the D*(2640) with our pole for the case oi S = 1 and not S = 
is going to be clear in the next section when the Dn channel is considered. In the case of 
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Squared amplitude for 1=1/2 and S=2 



6.xl0'^ 

^4.xio''; ; 

2.xl0'^- 

o i I - 
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Figure 8: Squared amplitude for / = 1/2 and S = 2 including the convolution of the 
p-mass distribution. 

the D*(2640) the width quoted in the PDG is very small, F < 15MeV. According to our 
result after taking into account the p mass distribution, one obtains 3-4 MeV, see Fig. [71 
In the case of D;{2A60), the width quoted in the PDG is 43 ±4 MeV for the Df and 
37 ± 6 MeV for -D2^. Then, it is clearly not compatible with the width found here, which 
is zero, see Fig. [HI and we need to allow that the resonance decays to another channel. For 
this reason we are going to consider the Dn channel in the next section, which is below 
the threshold of pD*. A novelty in this work is that we have obtained a new resonance in 
the sector 1=1/2 and 5 = which does not appear in the PDG, see Table [H One should 
note that if one takes Eq. f[32|) instead of the Eq. f[3T|) with a cutoff of the order of the 
natural size Qmax = 1-1-2 GeV, the results are very similar to those of Table [H (difference 
around 1%), which is an indication of the stability of the results. 



Channel 


D*{2600) 


^^(2640) 


L>^(2460) 


pD* 


14.32 


14.04 


17.89 


toD* 


0.53 


1.40 


2.35 



Table 5: Modules of the couplings Qi in units of GeV for the poles in the S = 0, 1,2; 1=1/2 
sector with the channel pD* and uD*. 
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5 The TiD decay mode 

5.1 Evaluation of the TrD-box diagram 

In the previous section we have obtained the positions of the poles and obtained a small 
width for the states taking into account the finite width of p. Here we consider a dominant 
decay mode into the nD channel in order to give account of a finite width. Our starting 
point is the set of diagrams of Fig. [91 One needs the pyrvr and the D*ttD vertex which are 
provided within the same hidden gauge formalism [1], [2], used in Section 2, by means of 
the Lagrangian 

Cv^^ = -ig{V^[^,d^^) . (37) 
For the first diagram of Fig. [9] we have: 

X (^3 - 2qYef{P + h~ 2g)"e(2)(P + - 2qfef 



X 



g2 _ _|_ (J^^ _ _ ^2 _|_ 



(P — qY — m\) + ie {k^ — qY — + ie 

Using the approximation that all the polarization vectors are spatial, it is possible to write 
the above amplitude as 



■At^D) 4 f ^"^Q ic (1) (2) (3) (4) 

1 



X 



q2 _ ^2 _|_ ^}^^ _ qY _ ^2 _|_ 



[P — qY — m'j) + ie {k^ — qY — rn^ + ie ^ ^ 

This integral is logarithmically divergent and as in [4J we regularize it with a cutoff in the 
three-momentum of the order of the natural size, for which we take qmax = 1-2 GeV. The 
results do not change much if one takes a value around this. After performing the dq^ 
integral of Eq. (!39|) , one finds 

Vi^D) ^ 4 / (1) (2) (3) (4) (1) (2) (3) (4) (1) (2) (3) (4)\ 
y \ i ' i j j ' J 



1 



IStt^ JO \2ujJ \ki + 2uj J k2—uj — ujD + ie 

111 1 



X 



k1 — uj — Ud + ie k^ — 2u + ie k^ — 2u! + ie P^ — u — ojd + 



1 / 1 \^ 1 



X ( 1 -^fiP^) im) 

P^ + uj + ud {k^ + uj + UDj 2ujd 
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where 



/(pO)= 4 {-32A;°PVcui3((P°)2-2cu2-3cucui3 



UJ 



D) 



+ 2{klfP'^UD{{Py-^ujl-2uuD-u:l) 

+ (fc3)^(2tu3 - {P'^foJD + Scj^cj^ + 2ujujI + cu|)) 

+ 18cucu^(-(P°)' + uol) + cu'(-12 (P°)='cuD + 44cu|,)) 

- (A;°)2(16cu^ + muo^uoD + 74cu=^cu^ + cub((P°)' - uolf 

+ 32cu2a;,5(-(P°)2 + a;^)+a;(-6(P°)2a;^ + 6cu^))} (41) 



and UJ = yg^ + '^^j = y"?^ + ^^1)5 P° = fc? + A;2- ■'■^ -^l- we can see clearly 
the sources of the imaginary part in the cuts k2{kl) — — Ud = 0, k^{k^) — 2u! = 0, 
— uo — uoo = 0, which give rise to the decays D*^ tiD, p — tttt and pD*^ tiD 
respectively. As in [4], to take into account the p mass distribution, we make a simple 
approach. First we neglect the three-momenta of external particles {ki ~ 0, i =1, 2, 3, 
4) and approximate ~ ^3 ~ "^p ^-nd fc^ ~ ^4 ~ m^.. Then we find double poles of 
{l/{k^ — 2u + ie))"^ and (1/(^2 — ^ — + ■ These double poles are then removed by 
replacing 

/ 1 \ 2 1 1 

(42) 



\kl-2uj + ie) kl-2uj + iTp/4: k\-2uj ~iTp/A 

and so on, considering a finite width for p and D* . Here we set Tp = lAQ.2MeV and 
r^i. = 2.1 MeV (results are identical if we put Td* = OMeV). Once this is done, the 
other diagrams of Fig. [9] are calculated easily, which takes into account the decay into the 
irD channel. One must make a projection into a proper isospin and spin. For isospin, only 
/ = 1/2 is allowed, while for spin, S = and 2 are allowed. Decay into S* = 1 is forbidden 
since the parity of the pD* system for p in s wave is positive, which forces the ttD system 
to be in L = 0, 2. Since the n and D have no spin, the total angular momentum J is equal 
to L in this case. Therefore, only the O"*", 2"*" quantum numbers have this decay channel. 
This provides an explanation on why the D*(2640) does not have practically width and 
the 1^2(2460) has a bigger width. Finally, we obtain 

^(27r,/=l/2,5=0) ^ 20 F^''^) , 

^(27r,/=l/2,5=2) ^ gyi^D) ^ (^43 ^| 

where V^'^^^ is given by Eq. PUI) after removing the polarization vectors. As in we use 
a form factor for a off-shell tt in each vertex, which is 

/\2 _ ™2 

= > (44) 

with A = 1400, 1500 MeV. The real and imaginary parts of the potential for the contribu- 
tions that we have calculated are plotted in Figs. [TD] and [TTJ As we can see, the real part 
of the potential coming from the vrD-box diagram is much smaller than the real part of 
the potential coming from the contact plus exchange terms. 
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Figure 9: vrD-box diagrams 



Real part of the potential for 1=1/2 and S=0 Real part of the potential for 1=1/2 and S=2 
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Figure 10: Real part of the potential for / = 1/2; S* = 0; and / = 1/2; S* = 2;. 
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Imaginary part of the potential for 1= 1/2 and S=0 

— ;rDbox, A=1400MeV ; 

— ;rD box, A= 1500 MeV 



Imaginary part of the potential for 1=1/2 and S=2 
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Figure 11: Imaginary part of the potential for 1 = 1/2; 5 = 0; and I = 1/2; S = 2;. 
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Figure 12: Squared amplitude for S = and 5 = 2 including the convolution of the p-mass 
distribution and the irD-hox diagram. 

5.2 Results with 



In Fig. [12] we show the results when one introduces the amplitudes of the Tables [H El [3] 
and the given ones in Eqs. fH3|) and fHOj) in the Bethe-Salpeter equation Eq. fl29|) . As one 
can see, now the states for 5 = and 5 = 2 have larger width since they can decay to ttD 
also. We have found a width of 20 MeV for the Z)2(2460), which is about 50% of the total 
width quoted in the PDG [6]. One could also have the D*tt decay channel, which would be 
possible by means of a anomalous coupling but, as it was seen in [1], this leads to a smaller 
contribution than the other one. Also in the PDG the most important contribution comes 
from the ttD channel. For the case of the new state found the width obtained is 50 MeV 
with A = 1500 Me V^. In the next section we introduce new elements of phenomenology 
that help obtain a somewhat larger width. 
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5.3 Results with using a different form factor and the ex- 

perimental coupling constant gD*DTT 

In this section we would like to evaluate again the vrZ^-box diagram, but, taking into 
account the strong coupling Qd'O-k and a different form factor also for an off-shell pion. 
After the recent measurements by the CLEO Collaboration [22], it is known that the 
D* meson couples strongly to Dtx. The experimental value of this coupling turns out to 
be larger, almost by a factor two, than the value obtained from some of the theoretical 
predictions using different approaches of the QCD sum rule [231 l2l]- Within the hidden 
gauge formalism, the vertex D*^D^7i^ obtained from Eq. fl371) is 

{D*+ip)n-iq)\D\p + q)) = -2 g'^,^^q,e'' , (45) 

with g'ij*ij^ = m/)*/2/£) = 6.3, which is also smaller than the experimental value for this, 
9D*Dn — 8.95 ±0.15 ±0.95. In [25], the D*D7r form factor is evaluated using the QCD sum 
rule for a D or a TT off-shell. A parameterization for an off-shell pion in terms of a form 
factor 

= QD^D^e'''/^' with A = 1.2 Gel^ , (46) 

is taken, together with another form factor to account for off-shell D mesons, which we do 
not need here since we are concerned about the imaginary part of the box diagram where 
the D meson will be on-shell. In Eq. f l46l) is a four momentum square = q^"^ — (f . 
With these assumptions on the form factors, a value of go'Dn = '^q'd'D-k = 14.9 is obtained 
in [25], which would be in better agreement with experiment {gn'On = 7.45). 

For the first diagram of Fig. [9l we consider the g° component of the vr^ on-shell, 
hence go = ~ 769.4 Mel/ and (fc? - g°) ~ 6 Mel/, for pD* at threshold, in the 

approximation of momentum zero for external particles. This leads to {k^ — q^Y/K^ ~ 10~^, 
for values of A around 1 GeV . Then it is licit to use the three-momentum version of the 
form factor of the Eq. for an off-shell pion in each of the vertices, that is, we replace 
in Eq. (j^Oj) the factor g'^ by 

9lA&.f{e-'^'"'y . (47) 
with ^fpTTTT = T^p/'^f-K = 4.2 and g^'Dn — 8.95 MeV" (the experimental value), A ~ 1- 
1.2 GeV and q running in the integral. 

In Figs. [13] and [H] we show the real and imaginary parts of the potential using Eqs. (140|) 
and p7j) . As we can see, the real part of the vrD-box diagram is very small compared with 
those coming from the contact term plus vector exchange terms, and therefore we can 
ignore it, thus focusing our attention at the imaginary part. The imaginary part is now 
larger than that in Fig. [11], but is still comparable with the values quoted in the PDG 
[6] for the width. We show the |Tp in Fig. [15] for various cutoff parameters. The |Tp is 
similar to Fig. [T21 but now the width is slightly larger. In the case of A = 1 GeV, we obtain 
40 MeV for 5 = 2, very close to the value quoted by the PDG (43 ± 4 MeV), and 61 MeV 
for S = 0. For A = 1.2 GeV we obtain for 5* = 2 a width around 60 MeV. Therefore, the 
two form factors, Eq. (jUj) and Eq. fHTl) . provide reasonable values of the width, with a 
preference for the option using the experimental g^o*D-K value and A = 1 GeV in Eq. P71) . 
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Figure 13: Real part of the potential for / = 1/2; 5 = 0; and / = 1/2; 5 = 2. 




Figure 14: Imagiary part of the potential for 7 = 1/2; 5 = 0; and I — 1/2] S — 2. 
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Squared amplitude for 1=1/2 and S=0 Squared amplitude for 1=1/2 and S=2 
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Figure 15: Squared amplitude for S" = and S = 2 including the convolution of the p-mass 
distribution and the vrD-boxed diagram. 

6 Conclusions 

We have made a study of the {puj)D* interaction using the hidden gauge formalism. The 
interaction comes from contact terms plus p meson exchange in the t-channel. Of all spin 
and isospin allowed channels in s-wave, we found strong attraction, enough to bind the 
system, in J = 1/2, S = 0, J = 1/2,^ = 1 and / = 1/2,5 = 2. We also found that 
in the case of / = 1/2, S = 2 the interaction was more attractive, than in the other two 
cases, leading to a tensor state more bound than the scalar and the axial vector. The 
consideration of the p mass distribution gives a width to the three states, rather small in 
all cases. Consideration of the ttD decay channel, in an equivalent way to what was done 
in the case of the pp interaction going to tttt in [Ij, makes the widths larger and realistic. 
Yet, the smaller phase space available here makes this contribution relatively smaller than 
in the case of the pp interaction. We found that the tensor state obtained matches the 
properties of the tensor state 1^2(2460). We predict two more states with S* = and 
S* = 1, which are less bound than the tensor state. We find in the PDG the state D*(2640) 
without experimental spin and parity assigned, but we conjecture that this state should 
be the S = 1 state found by us because we could find a natural explanation for the small 
experimental width of this state. The other state nearly degenerate in energy with this 
one, but with spin S" = 0, would still be to be found. The results obtained here should 
stimulate the search for more D states in the region of 2600 MeV . 
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